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An experiment conducted in the Mediterranean Sea in April 1996 demonstrated that a time-reversal 
mirror (or phase conjugate array) can be implemented to spatially and temporally refocus an 
incident acoustic field back to its origin. The experiment utilized a vertical source-receiver array 
(SRA) spanning 77 m of a 125-m water column with 20 sources and receivers and a single 
source/receiver transponder (SRT) colocated in range with another vertical receive array (VRA) of 
46 elements spanning 90 m of a 145-m water column located 6.3 km from the SRA. Phase 
conjugation was implemented by transmitting a 50-ms pulse from the SRT to the SRA, digitizing 
the received signal and retransmitting the time reversed signals from all the sources of the SRA. The 
retransmitted signal then was received at the VRA. An assortment of runs was made to examine the 
structure of the focal point region and the temporal stability of the process. The phase conjugation 
process was extremely robust and stable, and the experimental results were consistent with theory. 
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PACS numbers: 43.10.Ln, 43.30.Vh, 43.30.Bp, 43.30.Hw, 43.30.Re [DLB] 


INTRODUCTION 

Phase conjugation is a process that has been first dem¬ 
onstrated in nonlinear optics 1 and more recently in ultrasonic 
laboratory acoustic experiments. 2 ' 3 Aspects of phase conju¬ 
gation as applied to underwater acoustics also have been ex¬ 
plored recently. 4-7 The Fourier conjugate of phase conjuga¬ 
tion is time reversal; implementation of such a process over 
a finite spatial aperture results in a “time-reversal mirror. 2,3 ” 
In this paper we describe an ocean acoustics experiment in 
which a time-reversal mirror was demonstrated. 

In nonlinear optics, phase conjugation is realized using 
high intensity radiation propagating in a nonlinear medium. 
Essentially, the incident radiation imparts its own time de¬ 
pendence on the dielectric properties of the medium. The 
incident radiation is then scattered from this time-varying 
dielectric medium. The resulting scattered field is a time re¬ 
versed replica of this incident field propagating in the oppo¬ 
site direction of the incident field. For example, the scattered 
field that results from an outgoing spherical wave is a spheri¬ 
cal wave converging to the original source point; when it 
passes through the origin it has the time reversed signature of 
the signal which was transmitted from that point at the origi¬ 
nating time. Clearly, this phenomenon can be thought of as a 
self-adaptive process, i.e., the process constructs a wavefront 
of the exact required curvature. (An alternative would be to 
use a concave spherical mirror with the precise radius of 
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curvature of the incident wavefront.) There is an assortment 
of nonlinear optical processes which can result in phase 
conjugation. 1 In acoustics, however, we need not use the 
propagation medium nonlinearities to produce a phase con¬ 
jugate field. 

Because the frequencies of interest in acoustics are or¬ 
ders of magnitude lower than in optics, phase conjugation 
can be accomplished using signal processing. As in the op¬ 
tical case, phase conjugation takes advantage of reciprocity 
which is a property of wave propagation in a static medium 
and is a consequence of the invariance of the linear lossless 
wave equation to time reversal. In the frequency domain, 
time reversal corresponds to conjugation invariance of the 
Helmholtz equation. The property of reciprocity allows one 
to retransmit a time reversed version of a multipath dispersed 
probe pulse back to its origin, arriving there time reversed, 
with the multipath structure having been undone. 8 ' 9 This pro¬ 
cess is equivalent to using the ocean as a matched filter since 
the probe pulse arrival has embedded in it the transfer func¬ 
tion of the medium. This process can be extended further by 
receiving and retransmitting the probe signal with a source- 
receiver array. Depending on the spatial extent of the array, 
the above process results in some degree of spatial focusing 
of the signal at the origin of the probe signal. 

A time-reversal mirror (TRM) can therefore be realized 
with a source-receiver array. The incident signal is received, 
time reversed, and transmitted from sources contiguous with 
the receiving hydrophones. The time reversal can be accom¬ 
plished in a straightforward way, for example, by using the 
rewind output of an analog tape recorder or by a simple 
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program that reverses a digitized segment of a received sig¬ 
nal. 

An acoustic TRM has already been demonstrated in an 
ultrasonic laboratory using an array of source/receiver trans¬ 
ducers (SRA). 3 The array length was 10 cm and a single 
4-MHz source was placed at a transverse distance of 5 cm 
together with another receive array. The single source trans¬ 
mitted a probe pulse which was received at the SRA; the 
received pulse was time reversed and retransmitted from the 
SRA and subsequently received at an array (with the same 
orientation as the SRA) near the single source. The results 
showed a 15-dB peak at the location of the source relative to 
sidelobes away from the probe source location. Note that this 
focal point was at a range one-half the size of the aperture. 

Phase conjugation (PC) or the implementation of a TRM 
in the ocean is relevant to recent trends in acoustic signal 
processing which have emphasized utilizing knowledge of 
the environment, e.g., matched feld processing (MFP). 10 
However, MFP requires accurate knowledge of the environ¬ 
ment throughout the propagation path, which of course is 
difficult or impossible to obtain. Phase conjugation is an en¬ 
vironmentally self-adaptive process which may therefore 
have significant applications to localization and communica¬ 
tions in complicated ocean environments. Although the “ef¬ 
fective’ ’ ocean environment must remain static over the turn 
around time of the PC process, ocean variability on time- 
scales shorter than the turn around time might be compen¬ 
sated for with feedback algorithms. However, an understand¬ 
ing of relevant ocean time scales vis a vis the stability of the 
PC process will be required. 

In this paper we describe an April 1996 experiment in 
which an acoustic TRM was demonstrated in the ocean. In 
this initial experiment, a focal range of about 100 times the 
SRA aperture was accomplished easily with a 445-Hz probe 
source, a water depth of the order of 100 m, and a focal 
range of about 6.3 km. Large focal distances are obtainable 
in the ocean because in a waveguide geometry, a SRA has 
images which increase its effective aperture. Hence, there is 
an advantage to having a waveguide geometry over a free- 
field environment as was used first in the ultrasonic labora¬ 
tory experiment. Measurements in this first low-frequency 
ocean experiment also suggest a temporal stability of the PC 
process which is longer than what was expected intuitively. 
Some quantitative results on this stability are presented. 

In the next section we review the relevant theoretical 
issues including some simulation results leaving the details 
to an Appendix and appropriate references. Section II de¬ 
scribes the experiment in which the TRM was demonstrated 
and Sec. Ill presents the results. 

I. BACKGROUND THEORY AND SIMULATION FOR 
THE TRM EXPERIMENT 

The theory of phase conjugation vis a vis ocean acous¬ 
tics already has been presented. 4-7 Here we briefly review 
salient issues using the basic geometry of the TRM experi¬ 
ment (shown schematically in Fig. 2). More detail on theory 
is given in Appendix A and some additional details concern¬ 
ing experimental equipment are given in Appendix B. 



FIG. 1. Location of phase conjugation experiment. A source/receiver array 
(SRA) was deployed in 125-m-deep water and cabled approximately 1 km 
back to a small island, Formica Grande, the northernmost island of For- 
miche di Grosseto (42° 34.6' N, 10° 52.9' E). A rf telemetered vertical 
receive array (VRA) was deployed in 145-m-deep water approximately 6.3 
km west of Formica. 

A. Experimental geometry 

The TRM experiment was performed off the west coast 
of Italy in April 1996 as indicated in Fig. 1. Figure 2 is a 
schematic of the experiment and indicates the types of envi¬ 
ronmental measurements that were made. The TRM was 
implemented by a 77-m source-receiver array (SRA) in 125- 
m-deep water which was hardwired to the Isola di Formica di 
Grosseto. The SRA consisted of 20 hydrophones with 20 
contiguously located slotted cylinder sources with a nominal 
resonance frequency of 445 Hz. The sources were operated 
at a mean nominal 165-dB source level. The received signals 
were digitized, time reversed, and after being converted back 
to analog form, retransmitted. A probe source (PS) of the 
same type used in the SRA was deployed from the NATO 
research vessel ALLIANCE. The probe source together with a 
hydrophone was also used in parts of the experiment as a 
transponder. The ALLIANCE also deployed a vertical 46 ele¬ 
ment receive array (VRA) spanning 90 m located 6.3 km 
from the SRA which radio telemetered all individual element 
data back to the ALLIANCE. Ideally, the PS should be in the 
vertical SRA-VRA plane to correspond perfectly to the 
simulations in the following section; for practical reasons the 
PS was placed a few tenths of a km out of this plane. This 
did not introduce a significant error because the bathymetry 



FIG. 2. Experimental setup of the phase conjugation experiment. 
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FIG. 3. Collection of sound-speed profiles from CTDs taken during the 
experiment. 


was flat in the vicinity of the VRA. Figure 3 is a collection of 
the sound speed profiles (SSP’s) obtained from the 
conductivity-temperature-depth probe (CTD) as an indica¬ 
tion of the variability over the duration of the experiment. 
The bottom sound-speed structure as determined from earlier 
experiments 11 is shown in Fig. 4(a). More details on the in¬ 
strumentation and processing are given in Appendix B. 


B. Overview of theory 

For simplicity in this subsection, we analytically sum¬ 
marize the basics of phase conjugation in a range- 
independent waveguide. The simulations and comparisons 
with experimental data presented later on in the paper will 
include range-dependent modeling. The source, PS, is lo¬ 
cated a horizontal distance R from the source/receive phase 
conjugate array, SRA. 

1. Harmonic excitation 

The acoustic field, G u {R\Zj ,Z ps ), at the j th receiver el¬ 
ement of the SRA from the point source PS in Fig. 2 is 
determined from the Helmholtz equation 12 [assuming a har¬ 
monic time dependence of exp ( — icot)] 

y 2 G u> (r;z,z ps ) + k 2 (z)G l0 (r;,z,z ps ) 

oj 2 

= -S(r-r ps )S(z~z ps ), k 2 (z)=~ i^-y, (1) 

where z is taken positive downward and r =(x,y). Letting r 
be the horizontal distance from the probe source, Eq. (1) has 
the far field, azimuthally symmetric normal mode solution 
for pressure given by 


G m (nz,z ps )= p{Zi J 87Tr)U2 exp( — i tt/4) 


xE 


U n {z ps )u n {z) 


k 


1/2 

n 


exp (ik„r). 


( 2 ) 


where u„ ,k n are the normal mode eigenfunctions and modal 
wave numbers obtained by solving the following eigenvalue 
problem with well-known boundary conditions: 12 



(a) Sound Speed (m/s) 



FIG. 4. Single frequency simulation of phase conjugation for the geometry 
of Fig. 2 for a probe source located at a depth of 40 m and a range of 6.3 
km. (a) Sound-speed profile. The density, p and attenuation, a (in dB/ 
wavelength) of the bottom two layers are also given, (b) Simulation for a 
20-element SRA. Note the sharp focus in depth, (c) Simulation for only the 
bottom 10 elements of the SRA. 


d 2 U 

-^T + [k 2 (z)-k 2 ,]u n (z) = 0. (3) 

The mode functions form a complete set (for simplicity we 
omit discussion of the continuous spectrum though a good 
approximation is to use a set of discrete mode functions ob¬ 
tained from a waveguide extended in depth and terminated 
by a pressure release or rigid boundary) 
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2 

all modes 


u n (z)u n (z s ) 

p(z 5 ) 


8{z-z s ), 


and satisfy the orthonormality condition 

f“ u m (z)u n (z ) _ 

I / _ \ 

JO p(z) 


(4) 

(5) 


where S nm is the Kronecker delta symbol. 

The received field at the source/receiver array (SRA) at 
range R from PS with source/receive elements at depths Zj , 
is G JR\z,j ,z. ps ). The phase conjugation process consists of 
exciting the SRA sources by the complex conjugate of the 
received field, G*(R',Zj )• The resulting acoustic field trans¬ 
mitted from the J sources satisfies the wave equation, 

j 

V 2 Pp C ( r ’Z) + k 2 (z)P pc (r,z)='E S(z-Zj)Gl(R;Zj,z ps ), 

j= i 

( 6 ) 


k n ’s are nearly constant over the interval of the contributing 
modes) with the result that P pc (r,z) K=! S(z~z ps ). Figure 4 is 
a simulation of the phase conjugation process using Eq. (7) 
for a probe source at 40-m depth and at a range of 6.3 km 
from a 20 element SRA as specified in Fig. 2 verifying the 
above discussion. Range-dependent bathymetry was used as 
the input to an adiabatic mode model 13 for the specific 
sound-speed profile taken from the ensemble of profiles in 
Fig. 3 and a bottom sound-speed structure shown in Fig. 4(a) 
which includes a low speed layer as has been ascertained 
experimentally. 11 Notice that the focusing in the vertical is 
indicative of the closure property of the modes. As a matter 
of fact, for an SRA with substantially fewer elements, we see 
that the focusing still is relatively good. For example. Fig. 
4(c) also shows a result for the bottom 10 elements of the 
SRA which are below the thermocline. 

2. Pulse excitation 


where the range r is with respect to the SRA. Using Green’s 
function theory, the solution of Eq. (6) is the volume integral 
of the product of the Green’s function as specified by Eq. (1) 
and the source term of Eq. (6). For a vertical line of discrete 
sources, the integral reduces to a sum over the source posi¬ 
tions, 

j 

Ppc(r,z;a>)='2j G<o( r ',Z,Zj)G*(R;Zj ,Z ps ), (7) 

7=1 


In this experiment a 50-ms pure-tone pulse with center 
frequency 445 Hz was used for the probe transmission. We 
can Fourier synthesize the above results to examine phase 
conjugation for pulse excitation. Here, in the context of this 
experiment, we remind the reader that phase conjugation in 
the frequency domain is equivalent to time reversal in the 
time domain. The j th element of the SRA receives the fol¬ 
lowing time-domain signal, given by Fourier synthesis of the 
solution of Eq. (1): 


where R is the horizontal distance of the SRA from PS and r 
is the horizontal distance from the SRA to a field point. 

Note that the magnitude squared of the right-hand side 
(rhs) of Eq. (7) is the ambiguity function of the Bartlett 
matched-field processor 10 (with an appropriate normalization 
factor) where the data are given by G (JJ (R;z,j ,Z ps ) and the 
replica field by G l0 (r\z,Zj )■ In effect, the process of phase 
conjugation is an implementation of matched-field process¬ 
ing where the ocean itself is used to construct the replica 
field. Or, alternatively, matched-field processing simulates 
the experimental implementation of phase conjugation in 
which a source/receive array is used. To demonstrate that 
P pc (r,z ) focuses at the position of the probe source, ( R,z ps ), 
we simply substitute Eq. (2) into Eq. (7) which specifies that 
we sum over all modes and array sources 


P P c( r ->Z -,0)) 

m n j 


U m (z)u m (Zj)u n (Zj)u n (z ps ) 
p(Zj) p{z ps ) \Jk m k„rR 


P(R,Zj\t)= J GjR;Zj,z ps )S(<o)e- ia>t d<o, (10) 

where S( oj) is the Fourier transform of the probe source 
pulse. This expression incorporates all waveguide effects, in¬ 
cluding time elongation due to multipath propagation. For 
convenience, take the time origin such that P{R,Zj\t) = 0 
outside the time interval (0,r). Then the time reversed signal 
that will be used to excite the /'th transmitting element of the 
SRA is P{R,Zj\T—t ) such that T> 2 t. This condition is 
imposed by causality; the signal has to be completely re¬ 
ceived before it can be time reversed. Then 

P(R, Zj ;T-t)= j GjR;zj,z ps )S(oj)e- ia>(T - t) dw 

= f [G*(R; Zj ,z ps )e i ‘ oT S*((o)]e~ ia>t do), 

(ID 


Xexp i(k m r — k„R). (8) 

For an array which substantially spans the water column and 
adequately samples most of the modes, we may approximate 
the sum of sources as an integral and invoke orthonormality 
as specified by Eq. (5). Then the sum over j selects out 
modes m = n and Eq. (8) becomes 


P pf (r,z;w)^2 

m 


p(z ps )k m \frR 


exp ik m (r — R). 


(9) 


The individual terms change sign rapidly with mode number. 
However, for the field at PS, r = R, the closure relation of 
Eq. (4) can be applied approximately (we assume that the 


where the sign of the integration variable, u>, has been re¬ 
versed and the conjugate symmetry of the frequency-domain 
Green’s function and probe pulse has been used. The quan¬ 
tity in brackets in Eq. (11) is the Fourier transform of the 
signal received by the /th SRA receiver element after time 
reversal and time delay. Hence there is an equivalence of 
time reversal and phase conjugation in their respective time 
and frequency domains. 

Noting that the bracketed quantity in Eq. (11) is the 
frequency-domain representation of the signal retransmitted 
by the / th element of the SRA, Fourier synthesis can be used 
to obtain the time-domain representation of the field pro¬ 
duced by the TRM. Using Eq. (7), 
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Ppc(r,z;t) = 


i f 

j= 1 J 


G ul {r,Z,Z i )G* 0 (R,z j \z ps )e 


X5*(w)^'“' dco. 


( 12 ) 


This expression can be used to show that the TRM produces 
focusing in time as well as in space. Focusing in time occurs 
because a form of matched filtering occurs. To understand 
this, examine the TRM field at the focus point [that is, take 
r = R, z = z ps in Eq. (12)]. Neglecting density gradients, reci¬ 
procity allows the interchange G w (R,z pi ,Zy) 
= G u) (R,Zj,Z ps )■ Then the time-domain equivalent of Eq. 
(12) is 


P pc (r,Z\t) 



/ 


Gt* +t"^P’Zi iZps) 


XG t ,(R,Zj,z ps )dt 


S(t"-t+T)dt", 


(13) 


where the time-domain representations of the Green’s func¬ 
tion and probe pulse are used. Note that the Green’s function 
is correlated with itself. This operation is matched filtering, 
with the filter matched to the impulse response for propaga¬ 
tion from the probe source to the jth SRA element. This 
operation gives focusing in the time domain, that is, it re¬ 
duces the time elongation due to multipath propagation. 8 The 
sum over array elements is a form of spatial matched filter¬ 
ing, analogous to that employed in the Bartlett matched-field 
processor. 10 In addition, this sum further improves temporal 
focusing as the temporal sidelobes of the matched filters for 
each channel tend to average to zero which also is analogous 
to broadband matched-field processing results. 14 Finally, 
note that the integral over t" in Eq. (13) is a convolution of 
each matched-filtered channel impulse response with the 
time-reversed and delayed probe pulse. As a consequence, 
this pulse is not matched filtered, for example, a linear FM 
up-sweep will appear as a down-sweep at the focus and will 
not be compressed. 

Figure 5(a) shows a simulation for a 50-ms rectangular 
pulse with center frequency 445 Hz for the same geometry 
used in Fig. 4(a) as received at the SRA and Fig. 5(b) shows 
the pulse as transmitted to a plane at a range of 6.3 km, the 
range of PS. Four sources were excluded from the simulation 
because these phones were not used in the experiment. Note 
the temporal focusing; that is, the 50-ms pulse disperses to 
about 75 ms at the SRA but the time reversed pulse received 
at the VRA is compressed (focused) to 50 ms as opposed to 
exhibiting even further time dispersion. On the other hand. 
Fig. 5(c) shows a pulse 500 m outbound of PS (i.e., the VRA 
is at the same location but PS is 500 m closer to the SRA). 
The pulse is not spatially focused and it is temporally more 
diffuse than the result for the focal spot. 


3. Properties of the focal region 

A detailed discussion of the spatial and temporal factors 
affecting the focus is given in Appendix A. The primary 
result is that the TRM focus is robust, provided the SRA 
adequately samples the field in the water column. First, the 
focus tends to depend primarily on the properties of the 





FIG. 5. Simulation of a 445-Hz, 50-ms transmitted pulse for the geometry in 
Fig. 2 for a probe source located at a depth of 40 m. (a) Pulse received on 
the SRA at range of 6.3 km from PS. There is a temporal dispersion of about 
75 ms and significant energy throughout the water column, (b) The focus of 
the time reversed pulse at the VRA. There is pulse compression back to the 
original transmitted 50-ms duration as well as spatial focusing in depth, (c) 
Vertical and temporal distibution for a pulse 500 m outbound of PS (the 
VRA is at the same location but PS is 500 m closer to the SRA). 


ocean near the focus and tends to be independent of (the 
possibly range-dependent) properties of the medium between 
the SRA and the focus. Temporal changes in the medium due 
to, for example, surface waves and internal waves degrade 
the focus, but this degradation will be tolerable if the average 
(or coherent) Green’s function is not severely reduced by 
these time variations. Generally, the shape of the focus is 
approximated by the field that a point source placed at the 
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a) c) 




b) d) 




FIG. 6. Experimental results for probe source PS and VRA at same range, (a) The pulse data received on the SRA for PS at depth of 40 m. (b) The data 
received on the VRA from the time reversed transmission of pulses shown in (a). The VRA is 40 m inbound from the focus as determined by DGPS. (c) The 
pulse data received on the SRA for PS at depth of 75 m. (d) The data received on the VRA from the time-reversed transmission of pulse shown in (c). The 
VRA is 40 m outbound from the focus as determined by DGPS. 


focus generates after nonpropagating modes are subtracted. 
Thus if absorption or scattering tends to eliminate high-order 
modes, the focus will be comprised of the remaining lower 
order modes and will be relatively broader. Very roughly, the 
vertical width of the focus will be equal to the water depth 
(or depth of the duct) divided by the number of contributing 
modes if the sound speed (in the duct) is not strongly depen¬ 
dent on depth. 

The TRM focus is also robust with respect to array 
shape 4 provided the shape does not change between the 
probe reception and time reversed transmission. This prop¬ 
erty makes it unnecessary to know the exact shape of the 
TRM array and offers a considerable advantage over conven¬ 
tional beamforming. 

II. EXPERIMENTAL DEMONSTRATION OF A TRM IN 
THE OCEAN 

An assortment of runs was made to examine the struc¬ 
ture of the focal point region and the temporal stability of the 
process. Here we will be reporting on three types of experi¬ 
ments (note that range refers to the distance from the SRA): 

(1) Demonstration of the time-reversal mirror (TRM) in 
the ocean. The probe source (PS) is moved from shorter 
range to a longer range past the VRA. At each PS range, it 
sends out a 445-Hz, 50-ms pulse on the even minute. The 
pulse is received at the SRA, time reversed and retransmitted 


five times (once every 10 s) starting at the odd minute. This 
signal is received at the VRA and data from all channels are 
recorded. Note that when PS is at the same range of the 
VRA, the data recorded at the VRA are a vertical slice of the 
focal range as indicated in the simulation for a harmonic 
source in Fig. 4. Figure 5(b) is a simulation of the expected 
results at that range. When PS is closer than the VRA, the 
VRA data correspond to a measurement beyond the focal 
range and vice versa when PS is beyond the VRA. 

(2) Stability of TRM. PS is at the VRA range which 
means that we are measuring the vertical profile of the focal 
region. A 50-ms, 445-Hz pulse is sent out once and the SRA 
retransmits the same time reversed signal every 10 s for an 
extended period. Here the goal is to determine how long a 
single probe signal remains a valid phase conjugate probe for 
the specific ocean environment and source location. These 
results are constrained by the limitations of the actual experi¬ 
ment. 

(3) Acoustic ping pong. The probe source with colo¬ 
cated receiver now acts as a tranponder. The SRA transmits 
a 50-ms water column filling signal to the transponder which 
is at a depth of 75 m. The transponder retransmits the re¬ 
ceived signal (no time reversal) to the SRA which then trans¬ 
mits the time-reversed signals from the full array. This com¬ 
mences an acoustic ping pong iteration between SRA and PS 
with PS acting as a transponder (SRT). 
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a) 


J114131342 (OUT-OF-FOCUS) 



TIME (sec) 


c ) 


J114130132 (OUT-OF-FOCUS) _ FIELD DATA 

°1 SD = 40 m 

SR = 6.02 km 
RR = 6.24 km 



0.0 0.1 0.2 0.3 

TIME (sec) 



FIG. 7. Out of focus data received on the VRA from the time-reversed transmission of pulses with PS at a depth of 40 m. (a) PS is outbound 600 m. (b) PS 
is outbound 200 m. (c) PS is inbound 200 m. (d) PS is inbound 500 m. 


A. Demonstration of TRM in the ocean 

The vertical receive array VRA was deployed at a range, 
determined by DGPS, of 6.24 km from the SRA and the 
probe source PS was deployed at two different depths, 40 m 
and 75 m. Figure 6 shows the pulse as received on the SRA 
and VRA for both source depths. The data at the SRA are a 
combination of signal and noise. A 233-ms window was 
digitized and time reversed for transmission to the VRA. 
When the VRA and PS have the same range (experimentally 
within 40 m by a DGPS measurement) to the SRA, we see 
the focusing as predicted in Sec. I for a probe source at 40 m 
depth and similar results for a probe source at 75-m depth. 
Clearly, we have implemented a time-reversal mirror focus¬ 
ing at the range and depth of the probe source. 

Figure 7 shows the result as we sweep through the focal 
point. Note that because of the way the experiment had to be 
performed, we are actually keeping the VRA fixed and 
changing the range of PS An alternative way to present the 
focusing effect which displays the sidelobes off the main 
peak is shown in Fig. 8. The solid line with circles is the 
nearest to the focal region. Here we see the sidelobes in the 
vertical becoming large as we move away from the focal 
region. 

B. Stability of TRM 

The variability of the sound-speed structure in the water 
column is indicated in Fig. 3 which contains a collection of 


sound speeds derived from CTD’s at different positions and 
times throughout the experiment. A thermistor chain placed 
at the position indicated in Fig. 2 reveals the varying tem¬ 
perature structure as shown in Fig. 9. In addition, there is 
information concerning wave heights from the waverider 
shown in Fig. 2. The time series of the rms waveheight is 



Normalized Energy Level (dB) 


FIG. 8. The energy over a 0.3-s window as a function of depth for various 
ranges from the focal region. The depth of the probe source was 40 m. + 
means VRA is outbound from the focus (PS). 
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FIG. 9. Thermistor chain data. The contours from the top down are 15, 14.6, FIG. 10. Surface waveheight measurements from the waverider. 

and 14.2 °C. 
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d) 



FIG. 11. Results on stability of the focal region, (a) Pulse arrival structure at VRA for probe source at 40 m depth averaged over 1 h. (b) Pulse arrival structure 
at VRA for probe source at 75 m depth averaged over 2 h. (c) Mean and standard deviation of energy in a 0.3-s window for 40-m probe source, (d) Mean and 
standard deviation of energy in a 0.3-s window for 75-m probe source. 
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a) 


Mean Conjugate Field Strength Simulation, Rough Ocean Surface 




Mean Conjugate Field Strength Simulation, Rough Ocean Surface 



FIG. 12. Simulation of vertical profile of the mean field at the focal range 
for different values of surface roughness, (a) Probe source at 40 m. (b) Probe 
source at 75 m. 



FIG. 13. Acoustic ping pong between a transponder at 75 m depth and a 
range of 6.24 km from the SRA. The waterfall plot shows the energy in a 
0.3-s window at the VRA (which is at the same range as the transponder) as 
a function of depth for each of the 15 round trips. There were two minutes 
between each round trip. 

17:07 and J114 18:47-19:47, respectively. Figure 11 shows 
the results of these runs. These plots indicate that the focus 
was considerably more stable for the deep probe source ver¬ 
sus the shallower probe source and that the focus is broader 
for the shallower probe source. 

Simulations using representative rms wave heights from 
Fig. 10 and the environment of the experiment with a normal 
mode rough surface mean field scattering theory 15 are shown 
in Fig. 12. The results indicate that surface scattering does 
not have a significant impact on the focal region for this 
particular environment. On the other hand, examination of 
the environmental data indicates that the probe source at the 
shallower depth was at the bottom of the thermocline where 
the water column variability was the greatest. As derived in 
Appendix A, we expect the focusing phenomenon to be most 
sensitive to the environment at the endpoints of the experi¬ 
mental geometry. The tentative conclusion is that the fluc¬ 
tuations in this case were caused by sound-speed fluctuations 
in the water column, but more analysis and finer sampled 
volume data are required. 


shown in Fig. 10. Although the time series of the environ¬ 
mental data does not have the temporal and spatial resolution 
for an exhaustive comparison of theory and data, the two 
stability data collection periods show qualitative agreement 
with a first order analysis of the nature of the fluctuations. 

Basically, as shown in Appendix A, theory predicts that 
the mean field dominates the focal region with fluctuations, 
being a diffuse phenomenon, becoming more apparent away 
from the focus. That is, if one considers the total field to be 
composed of a mean field and a fluctuating field, it is the 
mean field which has the coherence properties which pro¬ 
duce the focusing whereas the fluctuating field is a form of 
signal-generated noise. 

Two stability data collection periods for the probe 
source depths of 40 m and 75 m were made for 1 h and 2 h, 
respectively (the lengths of the runs were dictated by experi¬ 
mental circumstance). The Julian day and times of the sta¬ 
bility runs for SD=75 m and SD=40 m were J114 15:11 — 



FIG. 14. Sound-speed profiles. The solid line was the optimum profile from 
the inversion process. The dashed line is thermistor chain derived sound- 
speed closest in time to the data shown in Fig. 6(a),(b). 
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a) 



TIME (sec) 


b) 



FIG. 15. Backpropagation using data from the SRA for the probe source at 
a depth of 40 m. (a) From the inversion process, (b) From the profile mea¬ 
sured at the time of the experiment. 

C. Acoustic ping pong: iterative focusing 

The purpose of the acoustic ping pong experiment is to 
demonstrate that focusing can be iteratively improved. This 
has already been demonstrated and explained in earlier free- 
held multiscatterer, ultrasonic experiments. 16-18 Basically, 
since a TRM returns signals to their origin in proportion to 
their original relative strengths, repeating the process a sec¬ 
ond time will reduce the level of the focused held for the 
weaker signals versus the stronger signals, and so on. The 
theoretical explanation is in terms of eigenvalues and eigen¬ 
vectors of the time-reversal operator. Eventually, only the 
strongest signal (or that part of the held corresponding to the 
largest eigenvalue) is focused. 

In this experiment, ping pong was initiated and kept go¬ 
ing for 15 round trips. Figure 13 is a waterfall plot of the 
energy in a 0.3-s window of the pulses received on the VRA 
which was at the same range as the transponder. There are 2 
min between each round trip. 

These results show the increased focusing brought about 
by the iteration process. However, this single source result is 
not completely analogous to the free-space multiscatterer re¬ 
sults in Refs. 16-18. Rather, it depends on the particular 
TRM array-data eigenvector structure in the specihc wave¬ 
guide environment. A paper with a detailed explanation of 
this process is in preparation. 
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FIG. 16. Backpropagation using data from every other element of the SRA: 
10-element TRM. (a) Probe source at a depth of 40 m. (b) Probe source at a 
depth of 75 m. 

III. EXTRAPOLATING THE EXPERIMENTAL RESULTS 

We have demonstrated that a time reversal mirror 
(TRM) can be implemented in the ocean and that its perfor¬ 
mance is consistent with theory. In this section we use a 
combination of data and theory to gain some additional in¬ 
sight into the potential usefulness of this process. In particu¬ 
lar, we examine: 

(1) its potential as a tool for inversion; 

(2) whether a smaller aperture or few source/receiver ele¬ 
ments would still be effective for producing a TRM. 

Further, we use item 1 to help estimate the TRM perfor¬ 
mance of a smaller SRA. 

A. TRM applied to sound speed inversion 

Empirical orthogonal functions 19 (EOFs) about the mean 
of the profiles shown in Fig. 3 were constructed. It was then 
found through trial and error that the mean profile was suf¬ 
ficient to provide the optimum focusing using simulated back 
propagation from the SRA data. That is, the coefficient of the 
first term of an EOF expansion was negligible with respect to 
the expected accuracy of the sound-speed profiles. A plot of 
this result compared to “CTD 11” which was used in the 
simulations in Sec. I is shown in Fig. 14. This procedure is 
akin to matched field tomography ’ ’ except that more 
information is available because of the vertical array at the 
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b) 



FIG. 19. Energy strength of five-element SRA backpropagation to the VRA. 
(a) Probe source at 40 m. (b) Probe source at 75 m. 

focal distance. Also shown in Fig. 14 is the sound-speed 
profile taken closest in time to the experimental runs under 
discussion. 

Figure 15 shows backpropagation results initiated from 
SRA data using (a) the profile obtained from the inversion 
and (b) the profile taken at the time of the TRM experiment. 
Clearly, the single experimental profile does not represent a 
range-independent profile descriptive of the experimental 
acoustic results, whereas the profile derived from the inver¬ 
sion represents an adequate range-independent approxima¬ 
tion to the structure of the water column. These results are 
also meaningful in the context of mismatch in matched-field 
processing. The experimental results indicate that a matched 
field processor using the measured profile would not localize 
the source. 

B. Reduced and sparse aperture TRM 

A reduced aperture SRA would enhance the practicality 
of an ocean TRM. We have already shown through simula¬ 


tion of a harmonic source in Fig. 4(c) that we can expect the 
phase conjugation process to remain effective as a focusing 
procedure as aperture is reduced. We should be able to reli¬ 
ably estimate the focal properties of a TRM using data from 
a subset of source/receiver elements and simulations of the 
backpropagation using the effective soundspeed profile 
shown in Fig. 14 found from the inversion. Of course, direct 
measurement for the sparse arrays would best study this as¬ 
pect of the TRM, but such data were not taken in this experi¬ 
ment. 

Figure 16 shows the results of an adiabatic mode model 
backpropagation of time-reversed pulse data from every 
other element of the SRA. We see that for both PS depths the 
focal region remains prominent for the ten-element SRA. We 
also present some results for an assortment of five element 
arrays in Figs. 17 and 18. Figure 19 shows a prediction of the 
vertical profile of the energy strength of these results which 
use five element subsets of SRA elements. The results are 
extended in depth to show the fields near the boundaries. The 
key thing to notice is that there are some very small arrays 
which still produce significant concentration of sound in the 
desired focal region. This probe source depth-dependent re¬ 
sult has practical ramifications for active sonar system con¬ 
cepts in which one desires to minimize boundary reverbera¬ 
tion at the range of the target. These results are not 
conclusive for the 40-m probe source depth because that was 
the depth of more or less maximum variation of the sound- 
speed profile. Hence, the sound-speed inversion result used 
in the backpropagation calculation might be the cause of the 
poorer focusing of the shallow source. 

IV. CONCLUSIONS 

We have constructed a time-reversal mirror (TRM) in 
the ocean and hence demonstrated that phase conjugation 
(PC) is realizable in the ocean using a source-receive array 
and rather simple signal processing. The waveguide nature of 
the ocean enhances the focusing properties over a free-space 
environment because the boundaries in effect enlarge the 
TRM aperture through its images. The degree of focusing is 
in excellent agreement with theory. Furthermore, an effective 
TRM need not be a full water column array. We also have 
investigated the stability of the PC process vis a vis ocean 
fluctuations and measurements suggest a relatively long sta¬ 
bility of the PC process. Future studies will be aimed at the 
detailed relationship between ocean variability and the PC 
process and an investigation into the possibility of using PC 
for inverting for the ocean environment. In addition, it 
should be straightforward to experimentally confirm predic¬ 
tions of the focal size versus SRA aperture. 
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APPENDIX A: FACTORS AFFECTING THE FOCUS 

In interpreting the results of the 1996 phase conjugation 
experiment, a primary issue is degradation of phase- 
conjugate focusing. Such degrading influences can be di- 
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vided into static and dynamic categories, the former includ¬ 
ing propagation and array structure effects and the latter 
including effects due to the time-varying ocean surface and 
volume. The object of study is the held produced by a phase 
conjugate source-receiver array (SRA), which can be writ¬ 
ten in the general form 

j 

Ppc(r,z;o)) = ^j G 2 (r,r J )Gf(r J ,r s ). (Al) 

j= i 

In Eq. (Al), P pc (r,z',a >) is the held produced at the held 
point, r =(r,z), by the phase-conjugate array with probe 
source placed at r ps = (R,z ps ) ■ The sum is over the J ele¬ 
ments of the SRA whose position vectors are denoted r ; - 
= (0 ,Zj). Following the convention used in the main text, 
horizontal ranges are measured from the SRA. Propagation 
from the probe source to the array is described by the 
Green’s function Gi(r„ ,r ps ), while propagation from the ar¬ 
ray to the held point is described by G 2 (r,r„). The sub¬ 
scripts 1 and 2 allow for the possibility that time variation of 
the ocean might cause changes in the Green’s function be¬ 
tween the probe and phase-conjugate transmission cycles. 
During either propagation cycle, the ocean is assumed to be 
“frozen” in the sense that it behaves as a time-invariant 
linear system. In this view, the Green’s function is the 
frequency-dependent system transfer function for acoustic 
propagation between any two points in the ocean. The fre¬ 
quency argument of the Green’s function used in the main 
text is suppressed here for convenience, but it becomes im¬ 
portant in treating pulsed transmissions. 

1. Phase conjugation in static environments 

The factors that control phase-conjugate focusing in 
static environments will be examined by considering a gen¬ 
eral nonuniform, nonadiabatic waveguide. The conditions for 
“ideal” phase-conjugate focusing in such a waveguide will 
be derived and this will implicitly identify the factors that 
degrade focusing. To simplify the discussion, only vertical 
phase-conjugate arrays will be considered. The main objec¬ 
tive is to generalize Eq. (9) of the main text to the range- 
dependent case, using the approach given by Siderius et al. 22 
in connection with the “guide source” concept. In this ap¬ 
proach, small regions near the probe source and SRA are 
assumed to be range independent, but the larger region be¬ 
tween is allowed to have arbitrary range dependence in 
bathymetry and sound speed. Losses are neglected and will 
be discussed later in qualitative terms. 

The Green’s function for the probe field near the probe 
source is approximated using range-independent normal 
modes. 


G w (r,z;R,z pi ) 


^ a„(z ps )u„(R,z) eikn(R) \ r _ R \ ^ 
n V k n (R)\r-R\ 


Similarly, the Green’s function for the probe field at the SRA 
is written in the form 


G M (0,z,;R,Zp S ) = 2 


b n(z p S )u n (0,Z /) 


VU0)R 


e * n (0)R 


(A3) 


The modal eigenfunctions in the vicinity of the probe source 
and SRA are denoted u n {R,z) and «„(0,z), respectively. The 
corresponding eigenvalues are k n (R) and k n ( 0). These 
Green’s functions do not bear the subscripts 1 and 2 intro¬ 
duced earlier because a time-invariant environment is under 
consideration. The subscript co is used here in the same sense 
as in the main text. The mode amplitudes for the near-source 
Green’s function are 


^n(z p s) 



W n (.P-iZ p s) i 


(A4) 


and the mode amplitudes for the Green’s function near the 
SRA are given by the linear transformation 


Mz) = 2 U mn a„(z). (A5) 

n 

For convenience, it is assumed that there are the same num¬ 
ber of modes near the source and near the array, so that U mn 
is a square matrix. Cases for which these numbers are similar 
but not equal can be treated by discarding high-order modes. 
The matrix U mn includes any mode coupling that is due to 
the range dependence of the ocean and is defined in such a 
way as to be independent of source depth. Furthermore, to 
the extent that absorption loss in the water column and sea¬ 
floor can be neglected, U mn is unitary. 

The field produced by the SRA is 

j 

P pe (r,z;w) = 2 GJr,z-,0,Zj)G* o (0,Zj',R,Z ps ). (A6) 

i= i 

The Green’s function for propagation from the y'th array 
element to the field point ( r,z ) can be expressed in terms of 
the Green’s function for propagation in the opposite direc¬ 
tion by using reciprocity: 

p{z) 

GJr,z;0,Zj)= GJ0,Zj;r,z). (A7) 

In terms of mode amplitudes. 


Gjr,z;0,zj) 


PW V c„(z)u n (0,Zj) oikn(m 

p(zj) n \fkjfijr 


(A8) 


where the mode amplitudes, c n (z), are 


G„(z) = 2 U mn a n (z)e ,k " {R){r R) 


(A9) 


Note that the mode amplitudes, c„(z), are essentially the 
same as the b n (z), but with the source range coordinate 
shifted by r — R. 

Equations (A3) and (A8) can be inserted in Eq. (A6) to 
obtain an expression for the phase-conjugate field in a range- 
dependent waveguide: 

p U,Z-,C0 )=-^ 2 Cm( z)AmlA * U > ) e /[A,„(())-*>)]« 

^Rr m,n ylk m (0)k*(0) 

(A10) 

where 


37 


J. Acoust. Soc. Am., Vol. 103, No. 1, January 1998 


Kuperman et al.: Phase conjugation in the ocean 


37 



A mn 


u m (0,Zj)u n (0,Zj ) 
P(Zj) 


(All) 


In the ideal case, the array spans the entire water column 
with elements having uniform spacing, d a , and the modal 
eigenfunctions have negligible amplitude in the bottom. In 
this case, the sum over array elements in Eq. (All) approxi¬ 
mates the orthogonality integral for modal eigenfunctions 
[Eq. (5)], and A,„„ d a can be taken equal to S m „ . This ideal 
can be approached quite closely in the environment of the 
1996 experiment. Using the environmental parameters de¬ 
fined in Fig. 4(a), and considering only the first 12 modes, an 
array with 36 elements with spacing d a = 3.33 m and with the 
shallowest element 4.44 m below the surface gives diagonal 
elements in A mn d a that are within 3% of unity and off- 
diagonal elements that are of order 0.03 or less. The first 
mode is an exception; it has a small diagonal element as it is 
trapped in the first sediment layer and not adequately 
sampled by the array. This is of no consequence, as this 
mode is very lossy and does not contribute to propagation. 
The element placement of the actual array gives smallest 
diagonal elements of about 0.5 with a few off-diagonal ele¬ 
ments as large as 0.3. 

Returning to the derivation of the conditions for ideal 
phase-conjugate focusing, take A mn = S mn /d a in Eq. (A10) to 
obtain 


P pc (r,z-,(o) 


p(z) 

d a \lR r 


2 Q mn a m (z)a*(z ps )e ,k >" ( - R)(r R) , 

m,n 

(A12) 


where 


Qmn 


y UlmUln -23[t„(0)]Jg 

/ *,( 0) 


(A13) 


Losses due to absorption and scattering are detrimental to 
phase-conjugate focusing, as they cause attenuation of 
higher-order modes, yielding a blurrier focus than would be 
possible with lower loss. Furthermore, this blurring will in¬ 
crease as the range between the source and the array in¬ 
creases owing to the strong range and mode number depen¬ 
dence of attenuation. Thus in defining the ideal case, losses 
are set to zero and the mode coupling matrix, U mn , is taken 
to be unitary. If the mode dependence of k/( 0) in Eq. (A13) 
is neglected, 

Qmn= kJ o)’ (A14> 

and the phase-conjugate field for an ideal array in a lossless 
environment can be approximated as 


P pc{r,Z\(o) = ^j 


n (R,z)u n (,R,z ps )e ik n (R)(r - R] 

%TTp(z ps )k n {0)d a 4Rr 


(A15) 


Apart from inessential factors, this expression is the same as 
Eq. (9) of the main text which was derived for the range- 
independent case. Even though Eq. (A 15) represents the 
ideal case, it illustrates properties that actual phase-conjugate 
arrays may possess, provided they are not too far from ideal. 


One such property is independence of the focus pattern on 
the distance between the probe source and the array (when 
absorption can be neglected and apart from the cylindrical 
spreading factor HyRr). Even more strikingly, the focus 
field is independent of the (possibly range-dependent) envi¬ 
ronment between the focus and the array (see examples pre¬ 
sented by Siderius et al .“). That is, the focus depends only 
on the local properties of the water column and sea floor and 
is not affected by bathymetry or range-dependent water col¬ 
umn properties in the region between the array and the focus, 
provided the latter do not change appreciably during the two 
propagation cycles. This means that, in the ideal case, phase 
conjugation is not affected by time-invariant forward scatter¬ 
ing due to bathymetry, fronts, etc. It also implies that, in 
simulations of phase-conjugate focusing, it is important to 
accurately model the ocean in the vicinity of the focus, but 
less accuracy is required for the more distant parts of the 
propagation path. One important reservation must be added 
at this point. The derivation above is essentially two dimen¬ 
sional in that cross-range spatial variation of the ocean is 
neglected. Static out-of-plane scattering will degrade phase- 
conjugate focusing if one-dimensional vertical arrays are 
used. Planar or volumetric arrays of sufficient aperture, on 
the other hand, will not suffer due to static out-of-plane scat¬ 
tering. 

The invariance seen in the ideal case is similar to that 
predicted for an ideal, closed phase-conjugate surface array 4 
which produces a strictly invariant focal field that resembles 
the original field of the probe source, except that the phase- 
conjugate field is a standing wave. In the present case, the 
probe source field (including only propagating modes) is 
given by Eq. (A2) which can be put in the form 


- iir /4 


G(r ,z;R,z ps ) = 


p(z ps )^&Tr\r-R\ 


u n (R,z)u n (R,z ps )e' k " (R) \ r R \ 

x 2 j . . 

" VUfl) 

(A 16) 

Apart from a difference in spreading loss and an overall 
phase difference, Eqs. (A15) and (A16) are quite similar. 
There is a slight term-by-term difference owing to differing 
factors involving modal eigenvalues, but the primary differ¬ 
ence is in the propagation phase factor. The source field 
propagates away from the source location while the phase- 
conjugate field propagates past the source location in the 
direction away from the array. 


2. Phase conjugation in time varying environments 

Time-dependent forward scattering due to surface and 
internal waves causes change in the propagation characteris¬ 
tics of the medium in the time interval between the probe and 
phase-conjugate transmission cycles with attendant degrada¬ 
tion of phase-conjugate focusing. 5 In discussing scattering 
from a general point of view, it is convenient to decompose 
the Green’s function into coherent and incoherent parts: 

G a (r,r') = G(r,r') + SG a (r,r'). (A17) 
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The subscript a takes on the values 1 and 2 for the probe and 
conjugate transmission cycles, respectively. The coherent, or 
mean. Green’s function, G(r,r') is not assigned a subscript 
because the random time variations are assumed to be sta¬ 
tionary in the statistical sense. It will be assumed that suffi¬ 
cient time has elapsed between the probe and conjugate 
transmission cycles that variations in the two Green’s func¬ 
tions are uncorrelated. 

(SG 2 (r d ,T c )SGf(r b ,r a )) = (SG 2 (r d ,r c )SG 1 (r b ,r a )) 

= 0. (A18) 

This condition was very likely satisfied in the 1996 experi¬ 
ment with respect to scattering by surface waves, which have 
correlation timescales on the order of seconds, while the time 
between transmission cycles was measured in minutes and 
hours. Internal waves have relatively long correlation time 
scales, but the longer transmission intervals (several minutes 
to a few hours) of the experiment were most likely sufficient 
to produce decorrelation of fluctuations in volume scattering. 

Combining Eqs. (Al), (A17) and (A18), the mean 
phase-conjugate field is 

j 

/V(r,z; W ) = 2 G(T,Tj)G*(Tj,T p3 ), (A19) 

i= i 

and the variance of the field is 
|P pf (r,z;co)| 2 -|P / , c (r,z;«)| 2 

J J 

= 2 2 [G(r,rj)G*(r,rj,)Kjj,(r ps ) 
j=i j r —i 

+ Gfr ; „,r ; )G*( r ps ,r y )K jj l (r) + K jj ,(r)K jjl (r ps )], 

(A20) 

where 

Kjj,(T) = (6G a (Tj,T)8G*(Tj' ,T)). (A21) 

The covariance, Kjji( r), is proportional to the correlation 
between the incoherent field at elements j and j ' of the array 
with a unit point source situated at r. In deriving Eq. (A20), 
free use was made of reciprocity (which allows interchange 
of the two arguments of the Green’s function) and stationar- 
ity (which means that SG , and SG 2 have identical statistics). 

Equations (A 19) and (A20) are general and include 
three-dimensional scattering (i.e., in-plane and out-of-plane 
scattering). They lead to two general conclusions regarding 
focusing in the 1996 experiment for those cases in which 
sufficient time elapsed between the two transmission cycles. 
First, the mean focus field, that is, the focus field averaged 
over many independent probe-conjugate-transmission cycles, 
is obtained by using the coherent Green’s function in place 
of the actual (random) Green’s function. Second, and most 
important, the field near the focus does not fluctuate appre¬ 
ciably, that is, it is well approximated by the mean focus 
field. This conclusion is supported by careful inspection of 
Eq. (A20), which shows that the variance of the phase- 
conjugate field is not localized near the focus, but is spread 
diffusely in range and depth. Thus near the focus, the mean 


field dominates, unless scattering is strong enough to dimin¬ 
ish the mean Green’s function to such a degree that focusing 
is essentially destroyed. 

To see that the field variance is unfocused, it is neces¬ 
sary to discuss each term in Eq. (A20) The first term can be 
viewed as being proportional to the intensity of a phase re¬ 
versed retransmission of the incoherent field produced by 
scattering of the probe transmission. This retransmission will 
be directed back toward the the scatterers responsible for the 
incoherent component of the probe field, and these are 
spread over the entire volume and surface of the ocean. Simi¬ 
larly, the second term is proportional to the intensity pro¬ 
duced at the source location by a coherent retransmission of 
the phase-reversed incoherent field produced from a fictitious 
source placed at the field point (reciprocity is being used in 
this interpretation). Again, this retransmission will be diffuse 
and will not peak as the field point approaches the source 
location. The last term in Eq. (A20) is more difficult to as¬ 
sess. It is a double sum over all array elements of the product 
of covariances due to sources placed at both the field point 
and probe source location. If scattering and propagation are 
very complicated in a spatial sense, these covariances will 
not be strongly dependent upon the source locations. That is, 
the incoherent field produced by these sources does not con¬ 
tain information on the source location. If this is the case, the 
covariances will be largely independent of r and r ps , and the 
last term of Eq. (A20) will not peak as r approaches r p ,. 

APPENDIX B: HARDWARE DESCRIPTION 

The phase conjugation (time reversal mirror) experiment 
was carried out in April 1996 off the northwest coast of Italy. 
As shown in Fig. 2, a source-receiver array (SRA) was de¬ 
ployed in 125-m-deep water and cabled approximately 1 km 
back to a small island, Formica di Grosseto (42° 34.6' N, 10° 
52.9' E). A if telemetered vertical receive array (VRA) was 
deployed in 145-m-deep water approximately 6.3 km west of 
Formica and used to measure the structure of the acoustic 
field across the water column. The R/V 
Alliance received this rf telemetered data stream and also 
deployed a source-receive transponder (SRT) (echo re¬ 
peater) which also was used as a probe source (PS). 

The vertical source array portion of the SRA consisted 
of 24 slotted cylinder sources spaced 3.33 m apart (total 
aperture 76.6 m). The sources have a resonance at approxi¬ 
mately 445 Hz and a 3 dB bandwidth of approximately 35 
Hz as shown in Fig. Bl. Thus the SRA sources were sepa¬ 
rated by approximately one wavelength at their center fre¬ 
quency. Each source was hardwired individually back to the 
transmit control system on Formiche di Grosseto via a mul¬ 
tiple twisted pair umbilical cable. The transmit control sys¬ 
tem synthesized the low-level analog signals for each source 
and these then were amplified prior to coupling onto the 
umbilical cable. Based on a nominal driving level of 100 
VRMS, the nominal source level of the transducers was 165 
dB re: /rPa. 

In addition to the vertical source array, the SRA in¬ 
cluded a colocated (i.e., physically strapped together) verti¬ 
cal receive array consisting of 48 hydrophones spaced half 
the separation of the source array transducers. The time se- 
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Driving Voltage 100 V rms 
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FIG. B1. Transmitting voltage response (TVR) versus frequency for one of 
the slotted cylinder source array transducers. 

ries from each array element was sampled at f s = 1.5 kHz 
using 24 bit A/D converters, multiplexed onto a single digital 
data stream, and cabled back to Formiche di Grosseto via a 
separate coaxial umbilical cable. 23 The shore-based digital 
data acquisition system archived the data stream and enabled 
capturing short segments of the array time series (from the 
24 hydrophones colocated with the source array transducers) 
for time reversal and retransmission by the transmit control 
system. Due to high-level contaminants observed in four of 
the time series, these channels were set to zero during the 
retransmission process. 

The rf telemetered vertical receive array (VRA) con¬ 
sisted of 64 hydrophones in a nested configuration over a 
90-m aperture. 24 A 46 element subset of these hydrophones 
with 2-m spacing was used to generate the results discussed 
in the main text. The time series from each array element 
was sampled at f s = 1.2 kHz, multiplexed onto a single digi¬ 
tal data stream, and sent via if telemetry to the R/V 
Alliance for both quick-look analysis and archival pur¬ 
poses. 

Last, the source-receive transponder (SRT) (echo re¬ 
peater) and probe source (PS) consisted of a slotted cylinder 
transducer identical to those used in the source array, and it 
was operated at the same nominal source level of 165 dB re: 
1 /jPa. When used as an echo repeater, the SRT included a 
separate receiving hydrophone to sample the acoustic field at 
the depth of the source. In this case, a short segment of the 
received time series containing a SRA transmission was cap¬ 
tured, amplified, and retransmitted (without time reversal). 
When used simply as a source, the SRT transmitted a 50-ms, 
445-Hz pulse which probed the multipath structure of the 
channel. In this case, the SRA received the temporally and 
spatially spread transmission, time reversed and amplified 
the 24 time series, and retransmitted them from the source 


array transducers. By allowing the R/V Alliance to tow the 
PS slowly through the range between the SRA and VRA, the 

focal region of the phase conjugation process could be stud¬ 
ied and these results are discussed in the main text. 
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